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We describe here the measurement of the flux of neutrinos created by the decay of solar 8 B by the 
Sudbury Neutrino Observatory (SNO). The neutrinos were detected via the charged current (CC) 
reaction on deuterium and by the elastic scattering (ES) of electrons. The CC reaction is sensitive 
exclusively to f e 's, while the ES reaction also has a small sensitivity to v^s and v T 's. The flux of 
He's from 8 B decay measured by the CC reaction rate is </) CC (^ e ) = 1-75 ± 0.07 (stat.)^Q (sys.) ± 
0.05 (theor.) X 10 6 cm" J s _1 . Assuming no flavor transformation, the flux inferred from the ES 
reaction rate is <j> ES (i/ x ) = 2.39 ± 0.34 (stat.)to'il (sys.) X 10 6 cm" J s _1 . Comparison of <f> cc (v e ) to 
the Super-Kamiokande Collaboration's precision value of ES (i^) yields a 3.3cr difference, assuming 
the systematic uncertainties are normally distributed, providing evidence that there is a non-electron 
flavor active neutrino component in the solar flux. The total flux of active S B neutrinos is thus 
determined to be 5.44±0.99 X 10 6 cm _2 s _1 , in close agreement with the predictions of solar models. 



1 Introduction 



Over thirtvvears of solar neutrino experi- 
ments 0&iaJ'i have demonstrated that the 
flux of neutrinos from all sources within the 
Sun is significantly smaller than predicted 
by models of the Sun's energy generating 
mechanisms EH. The deficit is not only uni- 
versally observed but has an energy depen- 
dence which makes it hard to attribute to 
astrophysical sources: the data are consis- 
tent with ajiegligible flux of neutrinos from 
solar 7 Be BliJ, though neutrinos from 8 B (a 
product of solar 7 Be reactions) are observed. 
A natural explanation for the observations 
is that neutrinos born as VeS change flavor 
on their way to the Earth, thus producing 
an apparent deficit in experiments detecting 
primarily z^ e 's. Neutrino oscillations — either 
in vacuum or matter — provide a mechanism 
both for the flavor change and the observed 
energy variations. 

While these deficit measurements argue 
strongly for neutrino flavor change through 
oscillation, a far more compelling demonstra- 
tion would not resort to model predictions at 
all but look for non-z/ e flavors coming from 
the Sun. The Sudbury Neutrino Observatory 



(SNO) was designed to do just that: pro- 
vide direct evidence of solar neutrino flavor 
change through the inclusive appearance of 
non-electron neutrino flavors from the Sun. 
We present here the first solar neutrino re- 
sults from SNO, which have also been de- 
scribed in an earlier publication c3. 



2 SNO Detector 

SNO is an imaging water Cerenkov detector, 
which uses heavy water (D2O) as both the 
interaction and detection medium t3. Fig- 
ure [l] shows a diagram of the detector. SNO 
is located ~ 2 km (6020 k.w.e.) under- 
ground in INCO Ltd.'s Creighton Mine, deep 
enough that the rate of cosmic ray muons 
passing through the entire active volume is 
just 3/hour. 

The 1000 tons of heavy water is con- 
tained in a 12 m diameter transparent acrylic 
vessel, and is surrounded by 2 ktons of light 
water shielding. The Cerenkov light pro- 
duced by neutrinos and radioactive back- 
grounds is detected by an array of 9500 
8 inch photomultiplier tubes (PMTs), sup- 
ported by a stainless steel geodesic sphere. 
Each PMT is surrounded by a light concen- 




Figure 1. Diagram of SNO Detector. 



trator, which increases the photocathode cov- 
erage to nearly ~ 55%. The front-end dis- 
criminator thresholds are set to fire on 1/4 of 
a photoelectron of charge. Outside the PMT 
support sphere is another 7 ktons of light wa- 
ter shielding. 

The detector is also equipped with a flex- 
ible calibration system, capable of placing 
sources almost everywhere in either the x — z 
or y — z plane. The sources that can be 
deployed include a diffuse multi-wavelength 
laser for measurements of optical parameters 
and PMT timing, a 16 N source which pro- 
vides a triggered sample of 6.13 MeV 7's, and 
a 8 Li source delivering /3's with an endpoint 
near 14 MeV. In addition, high energy energy 
(19.8 MeV) 7!s_are provided by a 3 H(p, 7) 4 He 
('pT') source 113 and neutrons by a Cf source. 
Some of the sources can also be deployed on 
vertical axes within the light water volume 
between the acrylic vessel and PMT support 
sphere. 

3 SNO Reactions 

SNO can provide direct evidence of solar neu- 
trino flavor change through comparisons of 
the interaction rates of three different pro- 



v e + d^p + p + e~ (CC) 
v x + d ^ p + n + v x (NC) 

The first reaction, the elastic scattering 
(ES) of electrons, has been used to detect 
solar neutrinos in other water Cerenkov ex- 
periments. It has the great advantage that 
the recoil electron direction is strongly corre- 
lated with the direction of the incident neu- 
trino (and hence the direction of the Sun). 
In addition, this reaction has sensitivity to 
all neutrino flavors. For z/ e 's, the elastic scat- 
tering reaction has both charged and neutral 
current components, making the cross section 
for v e 's 6.5 times larger than that for v^s or 

The deuterium in the heavy water makes 
the second process possible: an exclusively 
charged current (CC) reaction which (at solar 
energies) occurs only for v^s. In addition to 
providing exclusive sensitivity to z/ e 's, this re- 
action has the advantage that the recoil elec- 
tron energy is strongly correlated with the 
incident neutrino energy, and thus can pro- 
vide a good measurement of the 8 B energy 
spectrum. The CC reaction also has an an- 
gular correlation with the Sun which falls as 
(1 — 0.340 cos(# Q )) E-3, and has a much larger 
cross section (~ 10 times larger) than the ES 
reaction. 

The third reaction — also unique to heavy 
water — is a purely neutral current process. 
This has the obvious advantage that it is 
equally sensitive to all neutrino flavors, and 
thus provides a direct model-independent 
measurement of the total flux of neutrinos 
from the Sun. 

For both the ES and CC reactions, the 
recoil electrons are directly detected through 
their production of Cerenkov light. For the 
NC reaction, the neutrons are not seen di- 
rectly, but are detected in a multi-step pro- 
cess. When a neutrino liberates a neutron 
from a deuteron, the neutron wanders within 
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the D 2 and is eventually captured by an- 
other deuteron, releasing a 6.25 MeV 7 ray. 
The 7 Compton scatters an electron and it is 
this secondary particle which is detected. 

Although the data we present here were 
acquired with the acrylic vessel filled with 
pure D2O, the detector is now running with 
NaCl added to the heavy water. The addi- 
tion of the salt provides chlorine which has 
a larger capture cross section (and hence a 
higher detection efficiency) for the neutrons. 
The capture on chlorine also yields multiple 
7's instead of the single 7 from the pure D2O 
phase, which aids in the identification of neu- 
tron events. Eventually, discrete He 3 propor- 
tional counters will be added which will count 
neutrons exclusively. 

To determine whether neutrinos which 
start out as iVs in the solar core convert to 
another flavor before detection on Earth, we 
have two choices: comparison of the CC reac- 
tion rate to the NC reaction rate, or compar- 
ison of the CC rate to the ES rate. The for- 
mer has the advantage of high sensitivity — 
we compare the total flux to the u e flux 
and therefore expect to see a large differ- 
ence if the true neutrino flux agrees with 
standard solar models (which predict a total 
flux two to three times larger than previous 
measurements). In addition, uncertainties in 
the cross sections for the two processes will 
largely cancel. 

The second comparison has the advan- 
tage that both the CC and ES recoil electrons 
provide neutrino spectral information. The 
spectral information can ultimately be used 
to show that any excess in the ES reaction 
over the CC reaction is not caused by a differ- 
ence in the energy thresholds used to analyze 
the two reactions. The CC-ES comparison 
also has the advantage that the strong an- 
gular correlation with the Sun demonstrates 
that any excess seen is not due to some un- 
expected non-solar background. Lastly, the 
CC-ES comparison can be made with fairly 
high precision despite the small ES reaction 



cross section, because the Super-Kamiokande 
collaboration has already made a precision, 
high statistics measurement of the ES rate i. 
For the results presented here, only the CC- 
ES comparison will be described. 

4 Data Analysis 

The goal of the data analysis is the determi- 
nation of the relative sizes of the three signals 
(CC, ES, and neutrons) and ultimately the 
comparison of the rates. In the pure D2O de- 
tector configuration — the configuration with 
which these data were taken — we cannot sep- 
arate the signals on an event-by-event basis. 
Instead, we 'extract' the signals statistically 
by using the fact that they are distributed 
distinctly in the following three derived quan- 
tities: the kinetic energy of the recoil electron 
or capture 7 ray (T), the reconstructed radial 
position of the interaction (R 3 ), and the re- 
constructed direction of the event relative to 
the Sun (cos#©). 

Figure || shows these distributions for 
each of the signals. The top row of Figure || 
plots the different energy distributions for the 
three signals. We see in the figure that the 
strong correlation between the electron en- 
ergy and the incident neutrino energy for the 
CC interaction produces a spectrum which 
resembles the initial 8 B neutrino spectrum, 
while the recoil spectrum for the ES reaction 
is much softer. The NC reaction is — within 
the resolution smearing of the detector — 
essentially a (^-function, because the 7 pro- 
duced by the neutron capture on deuterium 
always has the same 6.25 MeV. 

The bottom row of Figure |^ shows the 
reconstructed direction distribution of the 
events. In the middle of that row we see the 
familiar peaking for the ES reaction, pointing 
toward the Sun. The ~ 1 — 1/3 cos 6q distri- 
bution of the CC reaction is also clear in the 
left hand side of the bottom row. Not sur- 
prisingly, the NC reaction shows no correla- 
tion with the solar direction — the 7 ray from 
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Figure 2. The energy (top), radial (middle), and di- 
rectional (bottom) distributions used to build pdfs to 
fit the SNO signal data. 



the captured neutron knows nothing about 
the incident neutrino. 

The distributions of reconstructed event 
positions is shown in the middle row of 
Figure pi These distributions are plotted 
as a function of R 3 , with R 3 = 1 occur- 
ring at the radius of the acrylic vessel (the 
edge of heavy water volume). We see here 
that the CC reactions — which occur only on 
deuterons — produce events distributed uni- 
formly within the heavy water, while the ES 
reaction (which occurs on any electron) pro- 
duces events distributed uniformly well be- 
yond the heavy water volume. The small 
leakage of events just outside the heavy water 
volume for the CC reaction is due to the res- 
olution tail of the reconstruction algorithm. 

The NC reaction, however, does not have 
a uniform distribution inside the heavy wa- 
ter like the CC reaction, but instead mono- 
tonically decreases from the central region to 
the edge of the acrylic vessel. The reason for 



this is that the capture cross section for neu- 
trons on deuterium is very small — the neu- 
tron wanders around long enough inside the 
D 2 that it may leak outside and be captured 
by hydrogen in either the acrylic vessel or the 
light water. Such hydrogen captures produce 
a much lower energy 7 ray (~ 2.2 MeV), be- 
low the analysis threshold. Therefore the ac- 
ceptance for events which are produced near 
the edge of the volume is reduced, because 
the probability of leakage there is correspond- 
ingly higher than for events produced near 
the center. 

One last point needs to be made regard- 
ing the distributions labelled 'NC in Fig- 
ure they represent equally well the detec- 
tor response to all neutrons, not just those 
produced by neutral current interactions, as 
long as the neutrons are produced uniformly 
in the detector. For example, neutrons pro- 
duced through photodisintegration by 7 rays 
emitted by U or Th chain daughters inside 
the D2O will have the same distributions of 
energy, radial position, and solar direction as 
those produced by solar neutrinos. In the 
analysis described here, no separation is done 
between these neutrons and those from the 
NC reaction. 

To determine the size of the three sig- 
nals, then, we use these nine distributions 
to create probability density functions (pdfs) 
and perform a generalized maximum likeli- 
hood fit to the same distributions in the data. 
There are, however, two principal prerequi- 
sites that must be satisfied before we can even 
begin this 'signal extraction' process. First, 
we need to process the data so that it is in 
a form we can use to do the fits. For exam- 
ple, we need to reconstruct the events to give 
us positions and directions that can be used 
to produce distributions, and we need to cal- 
ibrate the energy of each event. Even more 
importantly, we must be sure that the only 
signals present in the data are the three for 
which we are doing the fits — we have implic- 
itly assumed that the backgrounds are neg- 
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ligible. To accomplish this, we apply cuts 
to the data to eliminate backgrounds and we 
must ultimately demonstrate that any resid- 
ual backgrounds are negligible. 

The second signal extraction prerequisite 
is that the distributions used in the fitting 
process must be good representations of the 
detector's true response. In other words, we 
must build a model of the detector's response 
to the three signals which can be used to gen- 
erate the pdfs used in the fit. The model 
needs to reproduce the response at all places 
in the detector, for all neutrino directions, 
for all neutrino energies, and for all times. 
The last requirement is necessary because the 
detector's response changes over time due to 
things like failed PMT's or electronics chan- 
nels. 

The analysis we describe here therefore 
has three major components before the final 
fitting stage: the processing of data to re- 
move backgrounds, the building of a model 
to fit the data, and the demonstration that 
the residual backgrounds are small enough to 
use the signal-only model in the fits. 



4-.1 Data Processing 

We recorded the data set used in this analy- 
sis between November 2, 1999, and January 
15, 2000. Roughly 40% of the time during 
this period was taken up either by calibra- 
tion source runs or downtime caused by mine 
power outages. Of the remaining good data, 
we selected runs to analyze based on criteria 
which were 'blind' to the data itself — whether 
enough channels were live, whether calibra- 
tion sources were present, whether water as- 
says were being run, etc. After passing this 
run selection stage, no further run removal 
was allowed from the data set, and the fi- 
nal total livetime amounted to 241 days. Ap- 
proximately 30% of the data was put aside 
to serve as a blind test of statistical bias. As 
no significant differences were found between 
this sample and the other 70%, all subsequent 



discussion here refers to the full data set. 

During this time, the primary trigger 
threshold was set to fire on a ~ 100 ns coin- 
cidence of 18 PMT's each exceeding a chan- 
nel threshold of ~ 1/4 photoelectron. This 
trigger threshold corresponds to an energy 
of roughly 2 MeV. The trigger reaches 100% 
efficiency at 23 hit PMTs. In addition to 
the 100 ns coincidence, we ran simultane- 
ously with other triggers, such as a pre-scaled 
(1:1000) lower threshold (11 PMTs) trigger, 
a trigger on PMT pulse height sums, and a 
pulsed (random) trigger. 

The raw data set is far from the clean 
distributions shown in Figure ^. In par- 
ticular, the data is contaminated by instru- 
mental backgrounds arising primarily from 
PMT light emission ('flasher PMTs'), static 
discharges in the neck of the acrylic ves- 
sel, or electronic pickup. Although these 
instrumental backgrounds are very distinct 
from the neutrino signal, they occur at far 
higher rates: flasher events, for example, oc- 
cur roughly once each minute compared to 
the five to ten neutrino events we expect each 
day. We therefore developed a suite of low 
level cuts designed to remove instrumental 
backgrounds while losing a minimum of neu- 
trino events. These cuts were applied before 
any reconstruction of the data was done, and 
used only primitive information such as the 
PMT charge distributions, the raw and cali- 
brated time distributions, hits in veto tubes, 
and event-to-event time correlations. Fig- 
ure U shows the effects of the progressive ap- 
plication of these instrumental background 
cuts to the raw data set, illustrating the mul- 
tiple orders of magnitude reduction in the 
overall number of events. 

In any case in which such a large reduc- 
tion is obtained, the obvious question is what 
is the consequent reduction in good events — 
how much acceptance loss have we incurred 
by applying cuts which remove more than 
three orders of magnitude of the instrumen- 
tal backgrounds? To measure this loss, we 
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Figure 3. Effects of progressive application of instru- 
mental background cuts. 



Figure 5. Resolution in x for Li source data. 




Figure 4. Acceptance loss from low level instrumental 
background cuts, measured with calibration sources. 



used triggered calibration sources which pro- 
vided both samples of Cerenkov events and 
isotropic light events, and applied the same 
cuts to the source data as to the neutrino 
data. Figure [| shows the acceptance loss as 
a function of the number of hit PMTs, for 
16 N, 8 Li, and laser data. Although there is 
evidence of a bias at high energies (high num- 
ber of hit PMTs) , the overall scale of the loss 
is very small, ~ 0.5%. 

For events passing this first stage, we re- 
constructed the vertex position and direction 
of the particle using the calibrated times and 
positions of the hit PMTs. The reconstruc- 
tion algorithm begins with maximum likeli- 
hood fits using only PMT times, seeded by 
positions fixed to a grid throughout the de- 
tector volume. The best fit vertex from this 



grid-seeded procedure is then used as a seed 
for a second level of fitting which uses both 
the PMT times and their angular distribu- 
tion to simultaneously fit both the position 
and the direction of the event. The fitting 
process includes cuts on angular figures-of- 
merit which test both the quality of the fit 
and the hypothesis that the event is a sin- 
gle Cerenkov electron. Figure || shows the 
vertex resolution for electrons produced by 
the 8 Li source, which provides a localized 
(~ 5 cm) set of electrons with a broad spec- 
trum of energies. At 16 N energies, the vertex 
resolution is 16 cm and the angular resolution 
is 26.7°. 

For each event surviving the reconstruc- 
tion stage, we assigned an energy based on 
the hypothesis that the event was a single 
Cerenkov electron. While the number of hit 
PMTs by itself is directly related to the event 
energy, it must be corrected for the number 
of live channels online when the event was 
recorded, any change in the overall detector 
gain with time, and the optical effects of the 
intervening media between the Cerenkov pro- 
duction point (the event position) and the 
photon detection points (the hit PMTs). The 
optical corrections were calculated using in- 
situ measurements of the detector's optical 
properties (attenuation lengths, PMT angu- 
lar responses, etc.) and account for both the 
vertex position and the event direction. To 
minimize uncertainties associated with late 



hits (reflections, scattering, noise), the op- 
tical corrections use only prompt (in-time) 
photons by requiring the fitted time residu- 
als of the PMT hits to be within a narrow 
window around At = ns. 

Figure ^| shows the calibrated response of 
the detector to 16 N data. In Figure ^a we see 
the energy distributions for data taken with 
the source at the center and at R = 465 cm. 
The only corrections made here are for the 
number of live channels online, and therefore 
the shift in the mean of the two distributions 
is due to the different optical response at the 
two positions. Figure |^b demonstrates how 
the two distributions coincide once the opti- 
cal corrections are applied. 

With the event positions and directions 
fit and the energy calibrated, we passed the 
data through a final stage of cuts aimed 
at ensuring that the remaining events were 
consistent with Cerenkov light. We defined 
Cerenkov light with two orthogonal cuts: one 
which tests the narrowness of the timing dis- 
tribution, and one which tests the angular 
distribution of PMT hits. The former is done 
by cutting on the ratio of prompt (in-time) 
hits to the total number of hits in the event, 
and the latter by using the average angu- 
lar distance between hit PMTs in the event. 
Figure [?] shows three data sets distributed 
in these two variables: data tagged by the 
low level instrumental background cuts (tri- 
angles), Cerenkov data from the 16 N source 
(open circles), and neutrino data (closed cir- 
cles). The box used to define Cerenkov 
light is also shown, illustrating how both the 
source data and the neutrino data lie inside, 
while the instrumental backgrounds stay well 
outside. We required all data in the final sig- 
nal sample to lie within the box shown in Fig- 
ure @. 

The reconstruction quality cuts and the 
'Cerenkov box' cuts contribute to the overall 
acceptance loss, and we measured the scale 
of this loss along with the losses by the low 
level instrumental background cuts by us- 
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Table 1. Data processing steps. 



Analysis step 


No. of events 


Total event triggers 


355 


320 


964 


Neutrino data triggers 


143 


756 


178 


Mrit >30 


6 


372 


899 


Inst, bkgrnd cuts 


1 


842 


491 


Muon followers 


1 


809 


979 


Cerenkov box cuts 




923 


717 


Fiducial volume cut 




17 


884 


Threshold cut 




1 


169 


Total events 




1 169 



ing Cerenkov and laser calibration sources. 
The systematic uncertainties on these losses 
are associated with the calibration sources 
themselves (source reflectivity and shadow- 
ing), the low level electronic calibrations (for 
example, ADC pedestals) , and changes in the 
detector over time. Using the calibration 
source data, we find that the total loss for 
all cuts is 1.4±g;£%. 

For the final signal sample, we further 
restricted events to be within a fiducial vol- 
ume of 550 cm and have a kinetic energy 
T > 6.75 MeV. The fiducial volume restric- 
tion minimizes backgrounds associated with 
the acrylic vessel, light water, and PMTs, 
while the energy threshold reduces radioac- 
tive backgrounds and neutron events in the 
final signal sample. Table [l] summarizes the 
data processing, from the total number of 




events in the raw data set to the final sam- 
ple of 1169 events. The table also includes 
the effects of cuts which remove the spalla- 
tion products of cosmic ray muons. These 
cuts remove all events within 20 s of a par- 
ent muon. At this stage of the analysis, we 
have a data set which has been reconstructed 
and calibrated, and has had the majority of 
backgrounds removed. However, before fit- 
ting the resulting distributions, we still need 
to build a model of the detector's response, 
and demonstrate that the background re- 
moval has been successful enough that we can 
perform the fits using a signal-only model. 

4-2 Model Building 

The model of detector response we have used 
in this analysis takes as its inputs the physics 
of electron and 7 interactions in matter, the 
geometry of the detector, the behavior of 
the front-end data acquisition electronics and 
trigger, and — most importantly — the same 
measured optical parameters used in the en- 
ergy calibration described in the previous sec- 
tion. The model is a Monte Carlo simulation, 
which combines these inputs as well as the 
state of the detector as a function of time 
(the number of channels online, the overall 
energy scale determined by the 16 N source, 
etc.) to produce a predicted response func- 
tion for all event positions, directions, and 
energies. This response function is what we 




Figure 8. Deployment positions for 16 N source for a 
D2O scan and an H2O scan. 

use to create the pdfs which are ultimately 
used to fit the data. 

To ensure that the model is correct, we 
tested it against Cerenkov data representa- 
tive of the neutrinos we are trying to de- 
tect, for as many positions, directions, and 
energies as possible. The degree to which 
the model does not correctly reproduce the 
various measurements sets the scale for the 
systematic uncertainties on the predicted re- 
sponse function. 

Figure || depicts the positions inside the 
D2O and H2O for some of the 16 N scans. For 
the dependence on energy of the energy re- 
sponse, we compared 16 N data to pT data 
(6.13 MeV 7's to 19.8 MeV 7's). For the de- 
pendence on position and direction we com- 
pared different source positions and different 
sources — the Cf neutron source, for example, 
provides a very different event position distri- 
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Figure 10. Differences between predicted and mea- 
sured vertex resolution for the 16 N source as a func- 
tion of position in the detector. 



bution than the 16 N source does, and samples 
many more positions within the volume. We 
also tested the dependence on data rate by 
varying the rates for some of the calibration 
sources. Figure ^ summarizes the differences 
between the predicted energy response and 
the measured response for various sources as 
a function of time. The overall systematic 
uncertainty on the energy scale determined 
through these measurements is 1.4%. 

We performed the same kinds of tests for 
the prediction of the reconstruction accuracy, 
and Figure [To] compares the vertex resolution 
measured with the 16 N source at various po- 
sitions to the model prediction. There is a 
small systematic shift (~ 1 cm) between the 
two, but otherwise the model tracks the data 
well. The model prediction of the angular 
resolution agrees very well with the measure- 
ments made with the 16 N source, and has a 
negligible contribution to the overall system- 
atic uncertainty on the measured fluxes. 

Jj.,3 Backgrounds 

We are not quite ready yet to fit for the sig- 
nal amplitudes, because we must still demon- 
strate that the data is free enough from back- 
grounds to justify the use of a model which 
contains only signal distributions. There are 
three classes of background: the instrumen- 



tal backgrounds discussed in Section 4.1 , high 
energy 7 rays from the phototube support 



sphere and cavity walls, and low energy back- 
grounds from radioactivity both within and 
without the D2O volume. 

To measure the residual instrumental 
backgrounds, we used the 'Cerenkov box' 
cuts described earlier. The low level cuts 
aimed at reducing the instrumental back- 
grounds and the higher level Cerenkov box 
cuts are independent and orthogonal — and so 
the fraction of the instrumental backgrounds 
which lie inside the Cerenkov box is the same 
whether the instrumental background cuts 
identify the events or not. We therefore used 
the fraction of identified instrumental back- 
grounds which lie inside the box and multi- 
plied it by the number of events in the 'clean' 
data sample which lie outside. From this, we 
found the fraction of the clean data sample 
inside the Cerenkov box which may be due 
to instrumental backgrounds missed by the 
low level cuts. As a fraction of the final CC 
data sample, this is < 0.2%, small enough to 
ignore in the fit for signals. 

The determination of the high energy 
backgrounds was similar in principle, but 
here we had at our disposal calibration 
sources which provide triggered samples of 
high energy 7 rays. The 16 N source is nearly 
ideal for this measurement, as it acts as a 
triggered 'point source' of events which — 
with the exception of energy spectrum — look 
exactly like the background we are trying 
to measure. To use this source to mea- 
sure the backgrounds, we deployed it near 
where the backgrounds originate — out (and 
beyond) the detector's active volume. We 
then measured the ratio of the number of 
inward-going 7 events reconstructing just in- 
side the source position (the 'monitoring' 
box) to the number of events reconstruct- 
ing inside the 550 cm fiducial volume. With 
the number of events in the final data sample 
which reconstruct inside the same (but now 
spherically symmetric) monitoring box, we 
determined the number of background events 
which lie insde the 550 cm volume by mul- 



tiplying by the source-measured ratio. We 
explored systematic uncertainties by varying 
the monitoring box size, the deployment po- 
sition, and by using Monte Carlo simulation 
to explore the variation of the leakage with 
energy. The final limit on this source of back- 
ground measured in this way is < 0.8%. 

Low energy backgrounds originate from 
several sources: radioactivity in the heavy 
water, the acrylic, the light water, and the 
PMTs. Their typical energy is - 2 MeV, 
and our energy threshold of 6.75 MeV is high 
enough that the leakage can only come from 
the tail of the background energy spectrum. 
The small fiducial volume of 550 cm also 
greatly restricts the number of events from 
the PMTs, light water, and acrylic. To esti- 
mate the number of events from low energy 
backgrounds which leak above the signal en- 
ergy threshold or inside the fiducial volume, 
we used a combination of radioassays, encap- 
sulated U and Th calibration sources, and 
Monte Carlo simulation. Figure |ll| shows 
that the radioactivity in the heavy water — 
as determined by radioassays — is well below 
the original target values. At these levels, 
simulation shows that the tail of the back- 
grounds above our energy threshold is negli- 
gible. In the light water, assays also show 
that the the backgrounds are near or be- 
low target levels, but because these levels are 
still relatively high, we deployed calibration 
sources to measure the fraction that recon- 
struct within the 550 cm fiducial volume. Of 
the external sources of background, by far 
the largest is the radioactivity in the PMTs 
themselves. With calibration sources placed 
near the PMT sphere, we measured an upper 
limit on the leakage of the PMT radioactivity 
of < 0.2% of the final CC rate. 

5 Results and Implications 

We now have satisfied all the pre-requisites 
for doing a signal extraction: we have a clean 
data set in which the backgrounds are low 
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Figure 11. Water assay measurements of radioactiv- 
ity in D2O. 



enough to justify a signal-only fit, and a 
model which correctly predicts the response 
function of the detector as measured by cal- 
ibration sources. For the neutrino spectrum 
input to the model, we use an undistorted 8 B 
shape B0. 

The maximum likelihood fit to the 1169 
events in our sample gives us 975.4±39.7 CC 
events, 106.1±15.2 ES events, and 87.5±24.7 
neutron events, where the uncertainties given 
are statistical only. Figure |l2| shows the best 
fit to the distribution of event directions with 
respect to the Sun. The elastic scattering 
peak can clearly be seen, but with the avail- 
able statistics, only a hint of the slope of the 
CC electrons. 

To convert the CC and ES event numbers 
into fluxes, we need to correct for the accep- 
tance of the cuts, the energy threshold, and 
the fiducial volume restriction. We then need 
to normalize by the interaction cross sections 
and the number of deuterons and electrons 
inside the fiducial volume. For the CC cross 
section, we use the calculation of Butler et 
alt3, and do not include any radiative correc- 
tions. The radiative corrections may serve to 
increase the cross section by up to a few per- 
cent 113, and therefore decrease the measured 
value of the flux (and ultimately increase the 
significance of any difference between the CC 
and ES fluxes). We also include small correc- 
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Figure 12. Reconstructed direction distribution of electrons with respect to the Sun. 



tions due to the isotopic abundances of 17 
and 18 0, upon which CC reactions can also 
occur. Finally, we normalize by the overall 
livetime. 

Table || lists the systematic uncertainties 
on the flux measurements. The dominant un- 
certainties arise from our lack of knowledge 
of the true response function of the detec- 
tor. As described above, we characterize the 
scale of the uncertainties on the model by 
comparing the model predictions to measure- 
ments made with calibration sources, for ex- 
ample the 1.4% on the energy scale. To derive 
the uncertainties on the fluxes shown in Ta- 
ble ||, we varied the model predictions over 
the range of the uncertainties and repeated 
the analysis. In some cases this resulted in a 
larger uncertainty on the flux measurement — 
the 1.4% uncertainty on the energy scale be- 
comes a ~ 6% uncertainty on the flux derived 
from the CC rate, for example. 

In addition to the measurement of the 
systematic uncertainties, we have explored 
the systematic behavior of our results under 



Table 2. Systematic uncertainties on fluxes. 
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many different analysis approaches: for ex- 
ample, comparing different suites of low level 
cuts, reconstruction algorithms, Cerenkov 
box cuts, and choices of fiducial volume. We 
have also compared the results we get using 
the total number of hit PMTs as the measure 
of energy scale (thus changing the sensitivity 
to the knowledge of the late light distribu- 
tion) to the results from the energy calibra- 
tion described above. Lastly, we have per- 
formed fits using an analytical (as opposed 
to Monte Carlo) model of the detector re- 
sponse. In all cases, the results from these 
alternative approaches agree with the fluxes 
presented here to well within the systematic 
uncertainties quantified in Table ^. 

Converting the fit numbers to fluxes and 
including the systematic uncertainties listed 
in Table |2|, we find that the flux of neutrinos 
inferred from the ES reaction (assuming no 
flavor transformation) is 

^sno(^) = 2.39 ± 0.34(stat.) 

±g;12 (sys.) x 10 6 cm-V 1 

and the flux of 8 B f e 's measured by the CC 
reaction is 

0sno(^) - 1.75 ±0.07 (stat.)tSi? (sys.) 

±0.05 (theor.) x 10 6 cm^V 1 



where the theoretical uncertainty comes from 
the uncertainty in the CC cross section lij. 
The difference between these two numbers 
is 1.6cr, assuming that the systematic errors 
are distributed normally. The low signifi- 
cance of this result is driven mainly by the 
large statistical errors on the ES measure- 
ment. However, the Super-Kamiokande col- 
laboration has measured the flux with the ES 
reaction to high precision i, and finds 

0skK) = 2.32 ±0.03 (stat.)trlo? (sys-Xl) 
xlO 6 cm _2 s _1 . (2) 

The difference between SNO's measurement 
using the CC reaction (sensitive only to v^s) 
and Super-Kamiokande's measurement using 
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Figure 13. Summary of solar neutrino rate measure- 
ments from various experiments relative to BPB01 
standard solar model, including SNO and the derived 
8 B flux from the SNO and Super-Kamiokande rates. 



the ES reaction (sensitive to v^s and u T 's as 
well as v^s in the ratio of 1./6.5) is 3.3o\ This 
difference is therefore evidence of an active, 
non-z; e component to the solar 8 B neutrino 
flux. 

Figure |l3| summarizes the situation for 
all published solar neutrino experiments, in- 
cluding SNO. The points are plotted as ratios 
of the measured fluxes to the Standard Solar 
Model predictions of Bahcall, Pinnsoneault 
and Basu (BPB01 0), for the energy thresh- 
old used in each experiment. Here we can 
see the 3.3cr difference between the SNO and 
Super-Kamiokande measurements as well as 
the poor statistical accuracy of the SNO ES. 
Also plotted in Figure [l3| is the total flux of 
all 8 B neutrinos using the SNO CC measure- 
ment and the Super-Kamiokande ES mea- 
surement: 

4>{v x ) = 5.44 ± 0.99 x 10 6 cm-V 1 . 

We see in the figure that the agreement be- 
tween the measurement and the model pre- 
diction is very good. 



Figure |13| also shows that the differ- 
ences in the thresholds for the SNO and 
Super-Kamiokande measurements allows for 
the possibility that there is some spectral dis- 
tortion which could be causing the difference. 
Such a spectral distortion could occur if, for 
example, the oscillation were into a sterile 
neutrino. To look for such an effect, we can 
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Figure 14. Spectrum of electrons from v e interactions compared to prediction of the BPB01 standard solar 
model. 



first examine the spectrum of recoil electrons 
created by the CC interactions relative to the 
prediction for an undistorted 8 B spectrum. 
We derive such a spectrum by re-fitting the 
data energy bin-by-energy bin, without using 
the pdf for the CC energy spectrum. Fig- 
ure |lj shows the ratio of the spectrum de- 
rived this way to the standard solar model 
prediction. The dominant systematic uncer- 
tainties are indicated by the horizontal lines 
on the plot. Figure [li] shows that there is no 
large distortion in the expected spectrum. 

We can also eliminate the possibility of 
a spectral distortion leading to the difference 
in the SNO CC and Super-Kamiokande ES 
measurements by comparing the two mea- 
surements for the same neutrino energy. As 
described by Fogli et al E2l, this can be done 
by using different recoil energy thresholds for 
the SNO and Super-Kamiokande measure- 
ment. For these 'matched' thresholds (~ 
8.5 MeV for the Super-Kamiokande measure- 
ment compared to SNO's 6.75 MeV) we still 



get a difference of 3.1tr. This measurement is 
independent of both the standard solar model 
flux prediction and the predicted shape of the 
8 B spectrum — we need only know that the 
Sun produces 1/gS to see that there is a change 
to other active flavors. 

6 Future and Conclusions 

SNO's current and future data sets will pro- 
vide many more interesting measurements. 
We are now analyzing the pure D2O data 
in order to make a the measurement of the 
NC rate. The NC measurement should give 
us a confirmation of the CC-ES result, a 
higher precision measurement of the total 8 B 
flux, and a higher significance for the excess 
of non-z/g flavors (because we will be com- 
paring a 24 flux of 1.75 x 10 6 cm _2 s _1 to 
~ 5 x 10 6 cm~ 2 s -1 rather than to ~ 2.3 x 
10 6 cm~ 2 s _1 ). We are also analyzing the data 
in day and night bins, to determine whether 
any asymmetry is present (which would indi- 
cate that matter oscillations are the cause of 
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the flavor change as well as better restrict the 
allowed regions in the (tan 2 9, Am 2 ) plane. In 
addition, we are working on an analysis which 
includes the hep neutrinos. 

Beyond the pure D 2 data, we will also 
have the salt data set, which should provide 
us with an even better NC measurement, as 
well as new measurements of the other fluxes 
as well. Non-solar neutrino physics analy- 
ses are underway as well — looking at atmo- 
spheric neutrinos, anti-neutrinos (for which 
SNO has an exclusive coincidence tag), and 
supernova searches. 

SNO's first results, in combination with 
the Super-Kamiokande collaboration's mea- 
surements, provide direct evidence that solar 
neutrinos undergo flavor change on their way 
from the Sun to the Earth. They also show 
that the Standard Solar Model prediction of 
the 8 B flux is correct within the uncertainties 
of both the prediction and the measurement. 
We expect many more interesting measure- 
ments to come out of SNO for a long time to 
come. 
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